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On the basis of first-principles Go Wo calculations we study systematically how the electronic 
levels of a benzene molecule are renormalized by substrate polarization when physisorbed on different 
metallic and semiconducting surfaces. The polarization-induced reduction of the energy gap between 
occupied and unoccupied molecular levels is found to scale with the substrate density of states at the 
Fermi level (for metals) and substrate band gap (for semiconductors). These conclusions are further 
supported by GW calculations on simple lattice models. By expressing the electron self-energy in 
terms of the substrate's joint density of states we relate the level shift to the surface electronic 
structure thus providing a microscopic explanation of the trends in the Go Wo calculations. While 
image charge effects are not captured by semi-local and hybrid exchange-correlation functionals, we 
find that error cancellations lead to remarkably good agreement between the Go Wo and Kohn-Sham 
energies for the occupied orbitals of the adsorbed molecule. 

PACS numbers: 85.65.+h,71.10.-w,73.20.-r,31.70.Dk 



I. INTRODUCTION 

Solid-molecule interfaces are central to a num- 
ber of important areas of physics and chemistry 
including heterogeneous catalysis, electrochemistry, 
molecular- and organic electronics, and scanning tun- 
neling spectroscopyi£i£ii. Most of our current under- 
standing of level alignment at interfaces builds on effec- 
tive single-particle descriptions such as the Kohn-Sham 
scheme of density functional theory (DFT) 5 . Within such 
theories the energy levels of a molecule close to a surface 
are determined by hybridization, charge-transfer, and in- 
terface dipole fields - all properties of the static mean 
field potential defining the single-particle Hamiltonian. 
On the other hand, from photo-emission and electron 
transport measurements it is well known that the dy- 
namic polarizability of the molecule's local environment 
can have a large influence on the level positions^^^. 
Such polarization effects, which are induced by changes 
in the charge state of the molecule, are not captured by 
available single-particle descriptions. 

Many-body perturbation theory provides a system- 
atic method to obtain the true single-particle excita- 
tions [sometimes referred to as addition/removal ener- 
gies or quasiparticle (QP) energies] from the Green func- 
tion of the system. In the Go Wo approximation the 
electron self-energy is written as a product of the (non- 
interacting) Green function and a dynamically screened 
Coulomb interaction, E = zGoW o 11 ' 12 . It is instructive 
to compare this to the bare exchange self-energy given 
by H x = iGoV, where V is the unscreened Coulomb in- 
teraction. It is well known that the Hartree-Fock (HF) 
eigenvalues correspond to energy differences between the 
iV-particle groundstate and the unrelaxed N ± 1-particle 
Slater determinants (Koopmans' theorem). The effect 
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FIG. 1: (Color online) (a) Supercell used to represent benzene 
physisorbed on NaCl(OOl). (b) Reduction of a molecule's en- 
ergy gap when it approaches a polarizable surface, (c) Calcu- 
lated LDA, PBEO, and G W HOMO-LUMO gap of a ben- 
zene molecule lying flat at z = 4.5 A above different surfaces. 
Note that BaO(lll) is metallic due to surface states in the 
BaO band gap. 



of replacing V with the screened and frequency depen- 
dent Wo is two-fold: it introduces correlations into the 
many-body eigenstates, and it includes the response of 
the other electrons to the added electron/hole, i.e. re- 



2 



laxation effects. For a molecule at a surface, the latter 
effect is particularly important as it incorporates the at- 
tractive interaction between the added electron/hole and 
its induced image charge, into the QP spectrum. 

Recent experiments on molecular charge transport 
have renewed the interest for theoretical modeling 
of polarization-induced level renormalization. First- 
principles Go Wo calculations for a benzene molecule on 
graphite^ as well as CO on NaCl/Ge(001)— have demon- 
strated significant reductions of the molecular energy 
gap due to image charge effects. Model GW calcula- 
tions have been used to elucidate the qualitative features 
of the effect across different bonding regimes^. Classi- 
cal electrostatic models of various complexity have been 
developed to correct energy levels obtained from single- 
particle calculations ; 16 ' 17 ! 18 ' 19 . 

In this work, we present a systematic study of image 
charge-induced renormalization at a range of different 
surfaces taking both a classical and quantum many-body 
viewpoint. We have performed DFT calculations with 
local density approximation (LDA) and hybrid (PBEO) 
exchange-correlation functionals as well as Go Wo calcu- 
lations for a benzene molecule weakly physisorbed on 
the metals Li, Al, Ti, Rh, Pt, and the semiconduc- 
tors/insulators Ti0 2 , BaO, MgO, CaO, and NaCl. The 
results for the HOMO-LUMO gap of benzene are shown 
in Fig. [TJ While LDA and PBEO yields a substrate inde- 
pendent HOMO-LUMO gap, the Go Wo gaps are reduced 
from the gas phase value by an amount which depends on 
the polarizability of the surface. For all systems, we find 
that the dependence of the QP gap on the distance to 
the surface can be described by a classical image charge 
model. However, the model parameters are sensitive to 
the microscopic details of the system and this limits the 
usefulness of the classical model in pratice. By evaluating 
the Go Wo self-energy to second order we obtain a sim- 
ple analytic expression which relates the level shift to the 
substrate's joint density of states weighted by Coulomb 
interaction matrix elements. The model suggests that the 
HOMO-LUMO gap should scale with the substrate band 
gap (for semiconducting surfaces) and density of states 
at the Fermi level (for metallic surfaces). This trend is 
verified for the first-principles results and is further sup- 
ported by GW calculations for simple lattice models. Fi- 
nally, we analyze the deviation between the DFT and 
Go Wo results in more detail. We find that the occupied 
Kohn-Sham levels obtained with LDA (PBEO) are in very 
good agreement with the Go Wo results for benzene ad- 
sorbed on the metallic (semiconducting) surfaces, and we 
show that this is a result of significant error cancellation 
in the LDA/PBEO approximations. 

The paper is organized as follows. In Sec. |H]we outline 
the methodology used for the first-principles and model 
GW calculations. In Sec. IHII we investigate to what 
extent the first-principles Go Wo results can be explained 
by a classical image charge model. In Sec. HVI we derive a 
simple analytical expression for the polarization-induced 
level shift and show that it explains the main trends in 



both the first-principles as well as the model calculations. 
At the end of the section we analyze the description of 
occupied and unoccupied levels separately and discuss 
the effect of error cancellations in the DFT results. We 
conclude in Sec. [V] 



II. METHODS 
A. Ab-initio Go Wo calculations 

To model the solid-molecule interfaces we use a slab 
containing four atomic layers of the substrate in the ex- 
perimentally most stable phase, and a benzene molecule 
lying flat above the surface followed by 12A of vacuum. 
The benzene molecule is not relaxed on the surface but 
is fixed in its gas phase structure at a distance z from the 
surface. An example of a supercell is shown in Fig. []Ja) 
for the case of benzene on NaCl(OOl). The number of 
atoms included in the supercell per atomic layer is 9 for 
Al, Rh, Pt, Ti; 12 for Li and Ti0 2 ; and 16 for NaCl, MgO, 
CaO and BaO. This corresponds to distances between pe- 
riodically repeated benzene molecules in the range 8.1 to 
9.9 A. All DFT calculations have been performed with 
the PWSCF code^ 7 - which uses norm-conserving pseu- 
dopotentials^ 8 -. For exchange-correlation functionals we 
have used the local density approximation^ 9 , as well as the 
PBEO hybrid functional^ 1 . The Brillouin zone (BZ) 
was sampled on a 4x4x1 k-point mesh, and the wavefunc- 
tions were expanded with a cut-off energy of 40 Hartree. 

In the Go Wo (LDA) method one obtains the QP ener- 
gies from the linearized QP equation 



£ QP _ £ LDA 



LDA 



(e^)-v xc \^ A ) (1) 



where V'n DA an d e ™ DA are LDA eigenstates and eigenval- 
ues, and 



1 d(y£ DA |£ G w(6M DA ) 



ds 



1 -i 



-LDA 



(2) 



The self-energy, Egw, is evaluated non-selfconsistently 
from the single-particle Green function, i.e. Sqw — 
iG W , with G (z) = (z - H^)- 1 . It is customary 
to use the random phase approximation for the screened 
interaction, i.e. W = V{1 - VP)- 1 with P = -iG G . 

We have performed the Go Wo calculations with the 
Yambo code 3 - using the LDA wavefunctions and eigen- 
values from the PWSCF calculations as input. The plas- 
mon pole approximation has been applied with a fre- 
quency of 1 Hartree (the HOMO and LUMO energies of 
benzene change by less than 0.05 eV when the plasmon 
frequency is varied between 0.5 and 2.0 Hartrees). In the 
calculation of the self-energy we included a minimum of 
200 empty states. We have checked that calculations are 
converged with respect to slab thickness, lateral super- 
cell size, k-point mesh, all energy cut-offs, and that we 
reproduce the results previously reported in Ref. [H for 
benzene on graphite at z = 3.25 A. 
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B. Model GW calculations 

In addition to the first- principles Go Wo calculations, 
we have performed (self-consistent) GW calculations 
for two lattice models representing a metal-molecule 
and semiconductor-molecule interface, respectively. The 
model Hamiltonians contain three terms 



H = H so \ + i?mol + U, 



(3) 



describing the solid (metal or semiconductor), the 
molecule, and their mutual interaction, see Fig. O A 
metallic substrate is modeled by a semi-infinite tight- 
binding (TB) chain (we suppress the spin for notational 
simplicity) , 



Hmet = £(qci_i + cjljCi). 



(4) 



A semiconducting substrate is modeled by 
o 

H sc = ^2 ^2 £ af l ai+t(cl i c m - 1 +cl i _ 1 c m ), (5) 

a— c,v i— — oc 

where a — c, v refers to conduction and valence bands, 
respectively. 

The molecule is represented by its HOMO and LUMO 
levels, i.e. 



-Hmol = £,H fl H + 



(6) 



where e.g. fin — c jj-|- c K| + c hi° h I> is the number oper- 
ator of the HOMO level. 

Finally, the interaction between the molecule and the 
terminal site(s) of the substrate TB chain(s) is described 
by 



U = 



Un Nmoi 



for metals 



u Y,a( c to.a c vO,a + 4 CT c c0)CT )iV mo i , for semicond. 



where N mo \ = hn + is the number operator of the 
molecule. Note that since polarization of a semiconduc- 
tor occurs via transitions between valence and conduc- 
tion bands, only the interaction terms of the form given 
above contribute to the image charge effect (this will be- 
come clear in Sec. IIV A[) . 

We set Ep = corresponding to a half filled band for 
the metal. We choose and £l so that the molecule 
contains exactly two electrons (Ep in the middle of the 
HOMO-LUMO gap). We consider the limit of zero hy- 
bridization between the solid and molecule so that in- 
teraction between the solid and molecule occurs only via 
the non-local U . The model neglects interactions within 
the TB chain and between the molecule and interior TB 
sites (i < 0). These approximations are, however, not 
expected to influence the image charge physics described 
by the model in any qualitative way. 
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FIG. 2: (Color online) The lattice models representing a 
metal-molecule and semiconductor-molecule interface, respec- 
tively. We consider the weak coupling limit where no hy- 
bridization between the molecule and surface states occur. 
Thus the only interaction between the solid and molecule is 
via the non-local Coulomb interaction U. 



We obtain the Green function of the molecule from 

G(z) = l/(z-H inol -X GW [G}(£)) (7) 

where the Hartree potential due to U has been absorbed 
in ffmoi- The GW self-energy is calculated fully self- 
consistently using a recently developed GW scheme for 
quantum transport—. The renormalized molecular QP 
levels are obtained as peaks in the spectral function 
^(e) = -(l/7r)ImG^( £ ). 



III. CLASSICAL THEORY 

In this section we investigate to what extent the Go Wo 
results of Fig. [1] can be described by a classical image- 
charge model. 

The electrostatic energy of a point charge, q, located in 
vaccum at position (0, 0, z) above a polarizable medium 
filling the half-space z < zq, is given by (in atomic units) 



V 



qq 



4(z-z ) 



(8) 



The size of the image charge is q' = q(l — e)/(l+e), where 
e is the relative dielectric constant of the medium^!. In 
1973 Lang and Kohn showed that the energy of a clas- 
sical point charge above a quantum jellium surface fol- 
lows Eq. ([5]) with q' = —q (corresponding to e = oo as 
expected for a perfect metal), with the image plane, zq, 
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lying 0.5-0.9 A outside the surface depending on the elec- 
tron density^ 2 -. More recently, ab-initio Go Wo calcula- 
tions have found the same asymptotic form for the poten- 
tial felt by an electron outside a metallic surfac e 21 ' 23 ' 24 . 
From this it seems reasonable to conclude that the as- 
symptotic position of the electronic levels of a molecule 
outside a surface would also follow the image potential of 
Eq. ([5]). This is, however, only true for the unoccupied 
levels whereas the occupied levels experience a shift in 
the opposite direction, i.e. the shift is upward in energy 
as the molecule approaches the surface. This is because 
the occupied levels represent the negative of the energy 
cost of removing an electron from the molecule. Similarly 
it has been found that the image potential leads to band 
gap narrowing at semiconductor- metal interfaces . 25 ' 26 . 




FIG. 3: (Color online) Calculated Go Wo energy gap of ben- 
zene on NaCl, TiCh , and Ti surfaces (circles) as a function of 
the distance to the surface, and the best fit to the classical 
model Eq. (HJ (full lines). 



According to the classical image charge model all the 
molecular levels should experience the same shift (the 
sign of the shift being different for occupied and unoc- 
cupied levels). However, we have found that the best-fit 
values for zq and e mo dei obtained by fitting the HOMO 
and LUMO levels speparately, are in general different - 
most notably for the metallic surfaces. This observation, 
which is discussed in more detail in Sec. IIV C[ shows that 
the shape of molecular orbital also influences the size of 
the polarization-induced shift. 



TABLE I: Position of the effective image plane, zo, and di- 
electric constant, e m odel, obtained by fitting the ^-dependence 
of the HOMO-LUMO gap to Eq. © . Last row shows the ex- 
perimental optical dielectric constant of the bulk. The two 
values for the non-isotropic Ti02 refers to longitudinal and 
transverse polarization directions. 





zo (A) 


C mo del 


cxp,bulk 


NaCl(OOl) 


1.70 


1.15 


2.30 


MgO(OOl) 


1.20 


2.63 


2.95 


CaO(OOl) 


2.69 


1.56 


3.30 


BaO(OOl) 


2.74 


1.77 


3.83 


TiO 2 (001) 


1.79 


2.76 


8.43/6.84 


Al(lll) 


0.55 


oo 


oo 


Pt(lll) 


0.60 


oo 


oo 


Rh(lll) 


1.28 


oo 


oo 


Ti(001) 


1.66 


oo 


oo 


Li(001) 


1.72 


oo 


oo 


BaO(lll) 


2.01 


oo 


3.83 



Experimental data taken from Ref. |33 



To test whether the gap reductions obtained in the 
GoWo calculations can be described by the classical im- 
age charge model we have fitted Eq. ([8]) to the calculated 
HOMO-LUMO gap for z = 4.5, 5.5, 7.0, 9.0, oo A. In Fig. 
[2] we show the result of the fit for three systems (the 
fit is equally good for the other systems). The best-fit 
values for the effective image plane zq and the dielectric 
constant e mo dei are given in table Q] 

As can be seen e mo doi is generally smaller than the 
experimental optical dielectric constant of the bulk, 
e TC P,bulk - This is expected since the latter gives the long- 
range response of the bulk while e mo doi probes the local 
response at the surface. Part of the discrepancy between 
e™ p and e mo dci is clearly due to geometric effects. By 
taking the surface geometry into account, as done in Ref. 
[l8l . better estimates of e mo dci can be produced from e™ p . 
On the other hand, electronic effects due to the local 
atomic structure of the surface cannot be captured by a 
classical model. For example, the BaO(lll) surface is 
metallic due to surface states, and thus e mo dci ~ oo while 
e x P = 3.83. Similarly, impurities, defects, and surface 
roughness are expected to influence the local dielectric 
properties of the surface. 



IV. MICROSCOPIC THEORY 

In this section we first consider the GW self-energy for 
a molecule interacting with a surface to second-order in 
the electron-electron interaction. This leads to a sim- 
ple microscopic model for image charge renormalization 
which relates the shift of molecular levels to the electronic 
structure of the surface, and explains general trends of 
the first-principles and model GW calculations. In the 
last section we consider the HOMO and LUMO levels 
separately and explain how error cancellations in semi- 
local exchange-correlation functionals can explain the 
surprisingly good agreement found between LDA eigen- 
values and GW QP energies for the occupied levels of 
benzene on metallic surfaces. 



A. Second-order expansion 

In quantum many-body theory, the effect of substrate 
polarization on the energy levels of a molecule enters the 
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FIG. 4: (a) Feynman diagrams representing dynamic polar- 
ization of the substrate induced by an electron propagating 
in the molecule, (b) and (c): Generic shapes of the imag- 
inary and real parts of the self-energy of Eq. (|15f) for an 
occupied molecular level |a) interacting with a metallic and 
semi-conducting substrate assuming Vkk'.aa to be energy in- 
dependent. 



Green function via a self-energy operator. In general, the 
Go Wo self-energy can be written symbolically as 



E = = iG V(PVy 



(9) 



n=l 



where Go is the Green function of the non-interacting 
(Kohn-Sham) Hamiltonian, and P = —iGoGq is the po- 
larization bubble. The first-order term, E^, is sim- 
ply the static exchange potential while the remaining 
terms account for correlations and dynamic screening. In 
the following we consider the second-order term, E^ 2 ) = 
iGoVPV explicitly. This corresponds to approximating 
the reponse of the substrate by its non-interacting re- 
sponse, P. 

For sufficiently large surface-molecule separations (z > 
3.5A) we can neglect hybridization effects, and the non- 
interacting eigenstates of the combined system can be 
taken as the eigenstates of the isolated molecule and sur- 
face. We denote these eigenstates by {tpa} ('a' for adsor- 
bate) and {ipk}, respectively. To see how a given elec- 
tronic level, E a , is renormalized by polarization processes 
in the substrate we consider the (time-ordered) matrix 
element Eia'fw) = {->p a \T, (2 \uj)\ip a }, given by 

occ empty 

S ia = ^ ^ / iGo, aa (ui')V aa>kk < Pkk'(u'-Ul)V k >k,aadv'. 
k fc' 



(10) 

The Feynman diagram corresponding to Eaa is shown in 
Fig. 0Ja). The polarization and Coulomb matrices are 
given by 



'k/c' ,aa 



l 



i 



UJ - LO kk > +irj UJ + LU kk , - 17] 



(11) 



drdr' (12) 



where r\ is a positive infinitesimal and uj kk i = £k'—£k > 0. 
Using that G , Qa (w) = l/(a;-£ a +sgn(£ a - E F )ir])^, Eq. 
(fl"0l) reduces to 



7T J ' 



A(a/) 



uj — lu' + sgn(Ep — e a )ir] 
where we have defined the interaction strength 



duo' (13) 



occ empty 



A = TT^2 ^2 \Vkk>,aa\ 2 S(0J kk < - Sgn(e a - E F )(u - e a ))- 
k k' 



(14) 

Note that A is simply the joint density of states (JDOS) 
of the substrate, shifted by e a , and weighted by the 
Coulomb matrix elements. The physically relevant re- 
tarded self-energy is readily obtained from Eq. (flU)) 



<2hr = r r a(cj') 



TT UJ' — UJ 



-dui' — iA(u)). 



(15) 



where V denotes the Cauchy principal value. Now, the 
renormalized QP energy can be obtained from the equa- 
tion (neglecting off-diagonal terms) 



r QP 



(^|E( 2 )^( £ Q p )|Va}=0 (16) 



A graphical solution to the QP equation is illustrated in 
Fig. 0Jb,c) f° r the- case of an occupied molecular level 
E a < Ep interacting with a metal or semiconductor sur- 
face, respectively. 

From Eq. (fT4"|) it follows that the image charge 
effect does not broaden the molecular level because 
ImE( 2 )(£ a ) = 0. We also note that the level shift is in- 
dependent of the absolute value \s a — Ep\, and that the 
effect of changing the sign of £ a — Ep is to change the 
sign of the level shift. These properties are all in line 
with the classical theory. 

In the limit where V kk i, aa varies little with k and k' , A 
is simply proportional to the shifted JDOS [the "generic" 
cases illustrated in Figs. 0Jb,c)]. In this case the level 
shift is simply determined by the form of the JDOS. For 
a metal, the JDOS raises linearly at to = with a slope 
given by the metal's DOS at Ep. This suggests that 
the level shift should increase with the substrate DOS at 
the Fermi level. For a semiconductor, the JDOS raises 
smoothly at u> = E gapi suggesting that the level shift 
should decrease with E gap . In the following section we in- 
vestigate these relations for the model and first-principles 
calculations. We mention that the second order approx- 
imation discussed above may not always provide a good 
description of the full GW self-energy. However, as we 
will show in the next section, it explains qualitatively the 
trends in Go Wo calculations. 



B. Dependence of level shift on surface electronic 
structure 

In Fig. [5] we show the HOMO and LUMO levels of the 
lattice models calculated with the HF and GW approxi- 
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mations. In all plots we vary one parameter of the model 
while keeping the remaining parameters fixed2^. 
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FIG. 6: Calculated HOMO-LUMO gap of benzene at z = 4.5 
A(same numbers as in Fig. [TJ plotted as function of the LDA 
substrate band gap for the semiconductors (a), and the total 
DOS per volume evaluated at the Fermi level for the metals 
(b). Dashed lines have been added to guide the eye. (c) 
Average electron density in a plane lying z = 3.5 Aabove the 
clean surfaces. 



FIG. 5: HOMO and LUMO positions obtained from the 
simple lattice models for a metallic substrate (upper panel) 
and semiconducting substrate (lower panel). In all plots we 
vary one parameter while keeping the remaining parameters 
fixecOi. 

The upper panels refer to a metallic substrate and show 
the dependence on the levels on the interaction strength 
U and the intra-chain hopping parameter t. Note that 
the latter is inversely proportional to the projected den- 
sity of states (DOS) of the terminal site evaluated at Ep. 
The lower panels refer to a semiconducting substrate and 
show the dependence of the levels on U and the substrate 
gap, E gap . The HF eigenvalues are clearly independent of 
the non-local interaction between the molecule and sub- 
strate. This can be understood from Koopmans' theorem 
which states that the HF eigenvalues do not include the 
electronic relaxations of the substrate induced by the ex- 
tra electron/hole in the molecule. In contrast the GW 
levels vary in the way predicted by the simple model dis- 
cussed in the previous section: The polarization-induced 
reduction of the HOMO-LUMO gap is stronger for larger 
U as well as for larger substrate DOS at Ep for the metals 
and smaller substrate band gap for the semiconductors. 
A more detailed discussion of level renormalization based 
on the lattice model for metallic substrates, including the 
case of strong metal-molecule hybridization, can be found 
in Ref. QJ. 



In Fig. EKa,b) we plot the Go Wo gaps from Fig. [T] 
versus the LDA band gap and DOS at Ep for the semi- 
conducting and metallic substrates, respectively. For the 
semiconductors the reduction of the HOMO-LUMO gap 
clearly correlates with E gap . This indicates that the in- 
teraction strength, i.e. the matrix elements Vkk',aa of 
Eq. Q14p. do not differ too much from one surface to an- 
other. For the metals, the HOMO-LUMO gap seems to 
scale with the metal's DOS at Ep. However, we note that 
Li(001) and BaO(lll) deviate from the general trend fol- 
lowed by the other metals. This can be explained by 
the larger extend of the metallic wavefunctions of these 
systems into the vacuum region, which in turn leads to 
larger Vkk' ,aa matrix elements. Indeed, Fig. [6jc) shows 
the average electron density evaluated in a plane lying 
z = 3.5 Aabove the surface in the absence of the benzene 
molecule. The density outside the Li(OOl) and BaO(lll) 
surfaces is significantly larger than for the other surfaces 
which on the other hand have quite similar densities. 



C. DFT eigenvalues and error cancellation 

In Fig. [7] we plot the energies of the HOMO and 
LUMO levels of benzene at z = 4.5 A. For each sur- 
face, we have shifted the LDA, PBEO, and Go Wo levels 
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by the same amount so that the LDA HOMO is aligned 
with the HOMO in the gas phase. We note that the ef- 
fect of substrate polarization is very similar for the Go Wo 
HOMO and LUMO levels which are shifted up and down, 
respectively, by almost the same amount. This is indeed 
expected from the classical image charge model. Signif- 
icant deviations from this trend are, however, seen for 
Li(OOl) and BaO(lll). We ascribe this to the more ex- 
tended nature of the metallic states on these surfaces 
which reduce the validity of the point charge approxima- 
tion and can introduce differences between the Vkk' ,hh 
and Vkk' .ll matrix elements. 

Overall, the LDA and PBEO eigenvalues for the HOMO 
are in better agreement with the Go Wo QP energies than 
is the case for the LUMO. Moreover there is a general 
trend that the LDA eigenvalues come closer to the Go Wo 
energies as we move from the insulating to the metallic 
surfaces. In fact, the LDA HOMO level is almost on top 
of the Go Wo level on the metallic surfaces. This trend 
is clearly a result of sigificant error cancellation in the 
LDA. Indeed, it is well known that semilocal exchange- 
correlation functional overestimate (underestimate) oc- 
cupied (empty) molecular levels due to self-interaction 
effects. At the metallic surfaces this error is compensated 
by the missing image charge correction. PBEO gives bet- 
ter estimates for the free molecule where it opens up the 
LDA HOMO-LUMO gap due to partial removal of self- 
interaction errors. In this case, the cancellation between 
the missing image charge effect and the remaining self- 
interaction error results in very good agreement between 
PBEO and G W for the HOMO level on the semicon- 
ducting surfaces. 

The cancellation between self-interaction errors and 
missing polarization effects will always be present in 
hybrid- and semilocal approximations. However, the rel- 
ative size of the two contributions will in general depend 
on the shape of the molecule, its orientation with re- 
spect to the surface, the molecule-surface distance, and 
the type of substrate. 

V. CONCLUSIONS 

We have presented Go Wo calculations for a benzene 
molecule physisorbed on different metallic and semicon- 
ducting surfaces. Upon physisorption the molecule's 
HOMO-LUMO gap is reduced from its gas phase value 
due to dynamic polarization of the substrate. It was 
shown that a classical image charge model captures the 
qualitative features of the effect while the magnitude of 
the level shift is sensitive to the detailed atomic structure 
of the surface. In particular the presence of metallic mid- 
gap state at the surface of a semiconductor can have a 
large influence on the local response of the surface. Both 
local and hybrid exchange-correlation potentials fail to 
account for the polarization effects yielding Kohn-Sham 



eigenvalues of physisorbed benzene which are indepen- 
dent of the substrate. Nevertheless we found that a 




FIG. 7: LDA, PBEO, and G W energies for the HOMO and 
LUMO levels of benzene at z = 4.5 A above the surfaces. 
The very good agreement between LDA and Go Wo energies 
for the HOMO level at the metallic surfaces is due to error 
cancellation in the LDA approximation. 



cancellation between self-interaction errors and missing 
image charge effects in the LDA leads to a very good 
agreement between LDA and GoWo energies for the oc- 
cupied states of benzene on metallic surfaces. Similar 
conclusions were reached for the PBEO energies on semi- 
conducting substrates. Finally, we have derived a simple 
second-order approximation to the GW self-energy which 
expresses the polarization-induced shift of a molecular 
level in terms of the substrate's joint density of states 
weighted by Coulomb interaction matrix elements. This 
model was used to explain general trends in the first- 
principles results, namely the scaling of the benzene's 
HOMO-LUMO gap with the substrate density of states 
at Ep (for metals) and the substrate band gap (for semi- 
conductors). 

Our results clearly demonstrates the importance of 
non-local correlations for the electronic levels at solid- 
molecule interfaces. We expect this to have impor- 
tant implications for the theoretical modelling of electron 
transport in organic and single- molecule devices. 
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